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In the initial stage of relativistic heavy-ion collisions, strong magnetic fields appear due to the large ve¬ 
locity of the colliding charges. The evolution of these fields appears as a novel and intriguing feature in the 
fluid-dynamical description of heavy-ion collisions. In this work, we study analytically the one-dimensional, 
longitudinally boost-invariant motion of an ideal fluid in the presence of a transverse magnetic field. Interest¬ 
ingly, we find that, in the limit of ideal magnetohydrodynamics, i.e., for infinite conductivity, and irrespective of 
the strength of the initial magnetization, the decay of the fluid energy density e with proper time r is the same as 
for the time-honored “Bjorken flow” without magnetic field. Furthermore, when the magnetic field is assumed 
to decay ~ r _a , where a is an arbitrary number, two classes of analytic solutions can be found depending on 
whether a is larger or smaller than one. In summary, the analytic solutions presented here highlight that the 
Bjorken flow is far more general than formerly thought. These solutions can serve both to gain insight on the 
dynamics of heavy-ion collisions in the presence of strong magnetic fields and as testbeds for numerical codes. 

PACS numbers: 12.38.Mh,25.75.-q,24.85.+p,25.75.Nq 


I. INTRODUCTION 

Very intense magnetic fields of the order of B ~ 10 18 — 
10 19 G are produced orthogonal to the direction of motion 
in a typical non-central Au-Au collision at top RHIC energy 
(i.e., with a centre-of-momentum energy per nucleon pair of 
^/snn — 200 GeV). Recent studies show that the strength 
of the produced magnetic field grows approximately linearly 
with the centre-of-momentum energy of the colliding nucle¬ 
ons [ -3], It is now an experimentally well-established fact 
that in high-energy nucleus-nucleus collisions a very hot and 
dense phase of nuclear matter composed of quarks and glu¬ 
ons is formed. This hot and dense form of nuclear matter is 
also known as quark-gluon plasma (QGP). In the presence of 
a strong magnetic field as created in heavy-ion collisions, a 
charge current will be induced in the QGR leading to what is 
known as the “chiral magnetic effect” (CME) [ ]. At the same 
time, particles with the same charge but different chirality will 
also be separated, yielding what is called the “chiral separa¬ 
tion effect” (CSE). A density wave induced by these two ef¬ 
fects, called the “chiral magnetic wave” [4], is suggested to 
break the degeneracy between the elliptic flows of positive and 
negative pions [5], Moreover, it has also been found that there 
exists a deep connection between these effects and the Berry 
phase in condensed matter [6-8], Research on these topics is 
developing rapidly and a series of recent reviews and refer¬ 
ences can be found in Refs. [9- 1], 

The initial electric fields are also found to be quite large in 
event-by-event simulations of heavy-ion collisions [2], Such 
electric fields induce other novel effects such as the “chi¬ 
ral electric separation effect” (CESE) and the “chiral elec¬ 
tric wave” (CEW), which represent chiral currents and den¬ 
sity waves induced by the electric fields, respectively [12, 13], 
When the electric field is perpendicular to the magnetic field, 
just like for the Hall effect, a chiral current is expected, called 
the “chiral Hall separation effect” (CHSE), which might cause 
an asymmetric charge distribution in rapidity [ 4], 

Relativistic hydrodynamics has been proven to be quite suc¬ 
cessful in describing the experimentally measured azimuthal 


distribution of particle emission in non-central nucleus- 
nucleus collisions [15-20], It is then important to under¬ 
stand the effect of initial large magnetic fields on the fluid 
evolution. To this scope one needs a numerical code that 
solves the equations of (3+l)-dimensional relativistic mag¬ 
netohydrodynamics (MHD). There is consensus that due to 
very high velocities of the charges inside the colliding nu¬ 
clei (i.e., with Lorentz factors 7 ~ 100 for collisions at 
yGNN = 200 GeV), the magnetic fields decay very rapidly 
(i.e., decreasing by ~ 3 orders of magnitude within a timescale 
~ 1 fm for v /snn = 200 GeV Au-Au collisions) before the 
system reaches local thermal equilibrium and the fluid de¬ 
scription is applicable. However, the presence of a medium 
with finite electrical conductivity can substantially delay the 
decay of the magnetic field [21-23], 

Lattice-QCD simulations and theoretical models show that 
the QGP possesses a finite temperature-dependent electrical 
conductivity [24, 25], However, the interaction of the initial 
magnetic field with the QGP and its subsequent evolution is 
still an open issue and a topic of current research. An estimate 
of the relative importance of an external magnetic field on the 
fluid evolution can be obtained from the dimensionless quan¬ 
tity a = f? 2 /e, which represents the ratio of the magnetic- 
field energy density to the fluid energy density e. Clearly, val¬ 
ues of (7 > 1 indicate that one must consider the effect of the 
magnetic field on the fluid evolution. For a typical mid-central 
(i.e., with impact parameter b ~ 10 fm) Au-Au collision at 
top RHIC energy (-^/Snn = 200 GeV) the average magnetic 
field can be as high as ~ 10 ?n 2 ~ 10 19 G [.'., 2], where m n 
is the pion mass, which corresponds to an energy density of 
~ 5 GeV/fm 3 . Hydrodynamical model studies show that the 
initial energy density for such cases is ~ 10 GeV/fm 3 , thus 
implying a ~ 1 under these conditions. 

We note that the estimates made above are based on the as¬ 
sumption that the magnetic field (evaluated at time r = 0 fm) 
remains unchanged until the fluid starts expanding after reach¬ 
ing local thermal equilibrium at To ~ 0.5 fm. We also note 
that the estimate of a as given above is based on the event- 
averaged values for the initial magnetic field and energy den- 
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sity of the fluid. However, the situation can be very different. 
In fact, it is possible that the initial energy density distribution 
is very “lumpy”. Under these conditions, the produced mag¬ 
netic fields also show large variations and can be very large in 
some places where the corresponding fluid energy-density is 
small. In these cases, even for a quickly decaying initial mag¬ 
netic field we may locally have er > 1 up to the time when the 
hydrodynamical expansion starts. 

It is not the goal of this work to investigate the temporal 
evolution of the magnetic field produced in heavy-ion colli¬ 
sions. Rather, we concentrate here on the special case of one¬ 
dimensional, longitudinally boost-invariant fluid expansion a 
la Bjorken [26] under the influence of an external magnetic 
field which is transverse to the fluid velocity. In our analysis 
the evolution of the magnetic field is either regulated from the 
flux-freezing condition in ideal MHD or imposed in terms of 
a parameterized power law in proper time. The ultimate goal 
is that of finding analytic solutions for this flow that can be 
used both to gain insight in the dynamics of ultrarelativistic 
MHD flows as well as an effective test for more complex and 
realistic numerical codes (see also Ref. [28] for a work with 
similar intentions). 

The paper is organised as follows. Section II introduces our 
mathematical setup, while Sec. Ill presents the energy-density 
evolution when considering two representative prescriptions 
for the evolution of the magnetic field. A discussion of our 
main results is presented in Sec. IV, while a conclusive sum¬ 
mary is given in Sec. V. 

Following the predominant convention in relativistic hydro¬ 
dynamics of heavy-ion collisions, we use a timelike signature 
(+, and a system of units in which h = c = fc B = 1. 

Greek indices are taken to run from 0 to 3, Latin indices from 
1 to 3 and we adopt the standard convention for the summa¬ 
tion over repeated indices. Finally, we indicate three-vectors 
as bold face letter with an arrow and use bold letters without 
an arrow to denote four-vectors and tensors. 


II. MAT HE MATICAL SETUP 


We consider an ideal but magnetised relativistic fluid with 
an energy-momentum tensor given by 1 [29-32] 



where e,p, and a are the fluid energy density, pressure, and 
four-velocity, respectively. Since our considerations are re¬ 
stricted to special-relativistic flows, the metric tensor is that 
of flat spacetime, i.e., g 111 ' = rj^ v = diag(l, —1, —1, —1). 
Here B M F ua up is the magnetic field in the frame 

moving with the velocity up, e M " a/3 is the completely antisym- 


1 Note that expression (1) for the energy-momentum tensor is different from 

the one usually adopted in general-relativistic MHD (GRMHD) formula¬ 

tions and that we briefly review in App. A. 


metric four tensor, e 0123 = yj det \g\. The magnetic field four- 
vector B 1 ' is a spacelike vector with modulus B 1 ' B lt = —B 2 , 
and orthogonal to u i.e., = 0, where B = \B\ and 

B is the magnetic field three-vector in the frame moving with 
four-velocity 

As mentioned above, we are here interested in obtaining 
analytic solutions representing the MHD extension of one¬ 
dimensional Bjorken flow along the 3-direction with velocity 
uv = 7(l,0,0,u-) , where v z = z/t [26]. Hence, here¬ 
after we will assume the special case of a fluid flow in which 
the external magnetic field B is directed along the direction 
transverse to the fluid velocity v\ as remarked above, this rep¬ 
resents a rather good approximation of what happens in a typ¬ 
ical non-central Au-Au collision at top RHIC energy. This 
setup is also known as “transverse MHD”, since the magnetic 
field is contained in the transverse (x, y) plane [27], In ad¬ 
dition, since the fluid is expected to be ultrarelativistic, the 
rest-mass contributions to the equation of state (EOS) can be 
neglected and the pressure is simply proportional to the energy 
density, i.e., 

P = c 2 s e = - e , (2) 

where c s is the local sound speed which is assumed to be con¬ 
stant. The second equality in Eq. (2) refers to the case of an 
ultrarelativistic gas, or isotropic “radiation fluid” [ 33 ], where 
c s = 1/v/3 and which we will often consider in the remainder 
of this work. 

Rather than using a standard Cartesian coordinate system 
( t,x,y,z ), for longitudinally boost-invariant flow it is more 
convenient to adopt Milne coordinates, 

(U x, y, ?y) = ^\/f 2 - z 2 , x, y, ^ In j ) • ( 3 ) 

In these coordinates, the convective derivative is defined as 
u M dn = d T , while the expansion scalar takes the simple form 
0 = = t -1 . 

As customary, the projection of the energy-momentum con¬ 
servation equation = 0 along the fluid four-velocity, 

UyjdJT^ =0, (4) 


will express the conservation of energy. After some steps 
that can be found in App. A, we obtain the following energy- 
conservation equation 


dr 



e+p + B 2 

T 


= 0 . 


(5) 


Similarly, the projection of the energy-momentum conserva¬ 
tion equation onto the direction orthogonal to u, 

(Vuv - u lx u v )d ol T av = 0 . (6) 

leads to the momentum-conservation, or Euler, equation (see 
App. A), 

(e+p + B 2 ) d T u ll -(ry lv -Ui l u v )d v (p + =0. (7) 
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Note that for p = ip it reads 

|(H b2 ) =o ’ <s) 

thus showing that all thermodynamical variables depend only 
on r and are otherwise uniform in space. Considering instead 
p = x,y, Eq. (7) gives the MHD equivalent of the Euler equa¬ 
tion 


d-j-Ui 


1 


e + p+ B 2 


dA P + l -B 2 1=0. 


(9) 


With a uniform pressure and a magnetic field that depends 
only on r, the second term in Eq. (9) will vanish, thus imply¬ 
ing that if the velocities in the x- and //-directions are initially 
zero, they will remain so also at later times (i.e., d T Ui = 0). 


in. ENERGY-DENSITY EVOLUTION 

This section is dedicated to the discussion of two differ¬ 
ent cases for the evolution of the energy density depending 
on whether the magnetic field evolves according to the ideal- 
MHD limit (Sec. Ill A) or whether it follows an arbitrary 
power-law decay in proper time (Sec. IIIB). 


A. Ideal-MHD limit 

The solution of Eq. (5) requires the knowledge of the evo¬ 
lution of the magnetic field and hence of the induction equa¬ 
tion. In the limit of infinite electrical conductivity, i.e., in the 
ideal-MHD limit, the magnetic field obeys the frozen-flux (or 
Alfven) theorem and is thus simply advected with the fluid. In 
this case, setting B = \jB i B i , the induction equation takes 
the simple form 

B(t) = B 0 — , ( 10 ) 

Po 

where tq is taken to mark the beginning of the fluid expansion 
and po = p(7o), B 0 = B(tq), are the initial fluid rest-mass 
density and magnetic field, respectively. 

As written, Eq. (10) is of little use. In relativistic heavy-ion 
collisions, in fact, the net-baryon number is vanishingly small 
at mid-rapidity and the flux-freezing condition expressed by 
Eq. (10) needs to be modified to account for this. As we show 
in App. B, this is rather easy to do and yields, for an ultra- 
relativistic fluid with EOS (2), an evolution equation for the 
magnetic field, 

_ s -. /e\ 1/(1+c * } / e \ 3/4 

B(r) = B 0 - = B 0 - =B 0 - ,(11) 

so V e o / \ e o/ 

where s is the entropy density, so = s(ro), and the third equal¬ 
ity in Eq. (11) refers to the case c s = 1/v/3- Note that the sec¬ 
ond equality in Eq. (11) is the result of the first law of thermo¬ 
dynamics and reflects the relation between entropy and energy 
densities in an ultrarelativistic fluid [33], 


Using Eq. (11) in Eq. (5), we obtain 


B 2 0 


(£) 


2/(l+c® 


e+p + Bl (e/e 0 ) 2/(1+Cs) 


which can also be written in dimensionless form as 

e+p + er 0 e 2 /( 1+c ”) 


dr 


e+^e 2 /(i +c ») 


T 


= o, 


= 0 , 

( 12 ) 

(13) 


where e = e/eo, p = p/eo, and cr 0 = Bg/eo- 

Not surprisingly, in the limit of vanishing magnetic field, 
i.e., for /To —> 0, Eq. (13) takes the form of the Bjorken expan¬ 
sion, for which the energy density evolves according to 


d T e = = -(1 + c 2 )-, (14) 

T T 

where the second equality is written using the EOS (2). 

Using again Eq. (2), we find that both terms on the left-hand 
side of Eq. (13) contain a common factor that can removed: 
1 + (7oe( 1-c »)/( 1+c s)/(l + c 2 ). As a result, Eq. (13) can be 
written as 

d T e = -(1 + c 2 )- = , (15) 

T 6 T 

where the second equality refers to the more specific case of 
an EOS with c s = 1 / \/ 3. 

Thus, Eq. (15), which was derived within ideal MHD, co¬ 
incides with Eq. (14), which instead refers to Bjorken flow in 
the absence of an external magnetic field. After a more care¬ 
ful look, this result is not so surprising. In the ideal-MHD 
limit, in fact, the ratios B/s and thus B/e 1 ' /M fc 7 are con¬ 
served [cf. Eq. (11)] and although the total energy density will 
be larger in the presence of a magnetic field, the evolution of 
the fluid energy density will not be affected by the magnetic 
field which will be equally diluted as the fluid expansion takes 
place. This is essentially because the magnetic field has no ac¬ 
tive role in ideal MHD, but is simply passively advected in the 
expansion. Stated differently, a Bjorken flow is more general 
than formerly thought, as it applies not only to purely hydro- 
dynamical flows, but also to transverse MHD flows. To the 
best of our knowledge, this result, albeit natural, was not re¬ 
marked before in the literature. 

Equation (15) has an analytic solution of the form 

e (t) = cr _1/(1+c * ) = (y) ’ = (yj 1 , (16) 

where c is a constant that could be chosen, for instance, from 
the initial value of the energy density: eo = e( T o) = 1, and 
where the last equality again refers to c s = 1/ \J 'i. In the light 
of the remarks made above, it follows that Eq. (16) is also 
the solution for the energy density for Bjorken flow without 
magnetic field, cf. Eq. (14). 

Two final remarks: first, we note that combining Eqs. (11) 
and (16), it is easy to see that the evolution of the magnetic 
field in this case will be 

B(t) = ^P~ = e 1/(1+c * ) = £ 3/4 = ^ . (17) 

B o t 
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Second, our conclusion that the Bjorken flow is recovered in 
transverse MHD could have been reached also using entropy 
conservation and the Maxwell equations as long as the ther¬ 
modynamical relations are not affected by the presence of a 
magnetic field (zero magnetization vector), i.e., as long as 
de = Tds , where T is the temperature [34], In this case the 
derivation does not even require the specification of an EOS. 


B. Power-law decay 


Next, we explore cases where the external magnetic field 
does not vary according to the ideal-MHD flux-freezing con¬ 
dition (11), but has a different temporal evolution. Because we 
are in search of analytic solutions, we consider here a rather 
simple prescription and, in particular, one in which the mag¬ 
netic field follows a power-law decay in proper time, i.e., 

B(t) = B 0 (^) Q , (18) 

where a is a constant. Clearly, expression (18) is a simple 
ansatz but, as remarked in Eq. (17), it is sufficiently realis¬ 
tic to include the ideal-MHD case when a = 1. In addition, 
the range a > 1, i.e., of magnetic-field decay steeper than 
the ideal-MHD case, can be taken as a phenomenological de¬ 
scription of a resistive regime. Under these conditions, in fact, 
a finite electrical conductivity will lead to a more rapid decay 
of the magnetic field and, in turn, to a slower decay of the 
fluid energy density, which is “heated up” by the decaying 
field [cf. Eq. (20)]. 

Let us start by considering the equation of energy conser¬ 
vation (5), which for a general EOS of the form (2) and a 
magnetic-field evolution given by Eq. (18) yields 


dr 




o. (19) 


It is not difficult to find the analytic solution of this equation 
with initial condition So = 1, 



Once again, it is possible to see that in the limit of vanishing 
magnetization er 0 — > 0, Eq. (20) coincides with the solution 
(16) for Bjorken flow. Furthermore, for <to ^ 0 but a = 1, the 
solution (20) also coincides with Eq. (16), thus highlighting 
that Eq. (18) with a = 1 represents the evolution equation for 
a magnetic field in the ideal-MHD limit. 

Note that, for the second term in Eq. (20), the sign of the 
expression in brackets divided by 1+c 2 — 2a is always negative 
(remember that c 2 s < 1 by causality). Thus, for the case a > 1 
the second term is always positive. As a result, it always leads 
to a slower decay (and sometimes, as we will show below, 
even to an intermittent increase) of the fluid energy density 
than in the case a = 1. Viceversa, for a < 1 the second term 
in Eq. (20) is always negative, leading to a faster decay than 
in the case a = 1. 


Equation (20) seems to have a divergent behaviour at a = 
(1 + cl)/2, but this is only a first impression. We demonstrate 
in App. C that in the limit a —(1 + c 2 )/2, 

tm ( 7 q/ t ) 1+c " - = z z° y +c * i„ z!° z, 

a->(l+c2)/2 1+Cg — 2a V T / V T / 

( 21 ) 

Collecting things, the final solution of the energy- 
conservation equation (19) for a = (1 + c 2 s )/2 is 



Note that for r > r 0 the second term on the right-hand side 
of Eq. (22) is negative, hence increasing the decay of the fluid 
energy density with respect to the a = 1 case. Furthermore, 
while in the ideal-MHD case the solution e = 0 is obtained 
only asymptotically, for a = (1 + c 2 )/2 this extreme case is 
obtained after a finite time f = roe 2 ([( 1 _Cs ' )cro l. 

In summary, Eqs. (20) and (22) represent the solutions to 
Eq. (19); furthermore, since Eq. (20) comprises also the case 
a = 1, these equations provide a rather complete description 
of the full solution to the energy-conservation equation (19). 
As an example, we quote the solutions for c 2 = 1/3. For 
a/ 2/3, 



and for a = 2/3, 


(23) 



IV. DISCUSSION 

This section is devoted to a discussion of the various an¬ 
alytic solutions found in the previous section. For the sake 
of definiteness, we will always use the value c 2 = 1/3. We 
start by considering the ideal-MHD case, in which case the 
time evolution of the energy density and the magnetic field is 
given by Eqs. (16) and (17), respectively. These solutions are 
shown in Fig. 1 which reports the evolution of the normalised 
total energy density e + ^ao(B/B q) 2 for tq = 0.6fm. Dif¬ 
ferent lines refer to different values of the initial magnetiza¬ 
tion, ranging from <to = 0 (Bjorken flow without magnetic 
field; black solid line) up to cases with initial magnetization 
of (To = 1 (light-blue dotted line) and op = 10 (red dashed 
line). As already discussed in the previous section, the evolu¬ 
tion of the fluid energy density does not depend on op [cf. Eq. 
(16)] and scales like r -4 / 3 , while the magnetic energy density 
scales like t -2 . As a result, increasing op (as we do in Fig. 1) 
only adds energy density to the system, but does not alter the 
temporal evolution of the fluid energy density. 

Having considered the simple case a = 1, we next discuss 
the behaviour of the solutions when the magnetic field varies 
according to the more general power law (18). We have al¬ 
ready mentioned that a > 1 corresponds to the case when 
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r[fm] 

FIG. 1: Evolution of the normalised total energy density e + 
\oq(B/Bq) 2 . Different lines refer to different values of the ini¬ 
tial magnetization, ranging from ao = 0 (solid black line) up to 
cases with initial magnetization of ao = 1 (light-blue dotted line) 
and cro = 10 (red dashed line). Note that the fluid energy density 
decays like t~ 4 ^ 3 for all values of ao, i.e., as in traditional Bjorken 
flow. 


the magnetic field decays faster than in the ideal-MHD limit 
and could therefore be associated to a resistive regime. Con¬ 
versely, a magnetic-field evolution with a < 1 would corre¬ 
spond to a decay that is slower than in the ideal-MHD limit. 
As the case a = 1 is the limit of infinite conductivity, and 
thus of a maximal magnetic induction, it is at first sight hard 
to imagine how to produce a magnetic field that decays even 
slower than in the ideal-MHD case. However, in heavy-ion 
collisions the remnants of the colliding nucleons can give an 
additional contribution to the magnetic field, slowing down its 
decay [ 2 ]. Thus, considering also the case a < 1 is reasonable 
in this context. Within this range, a particularly interesting so¬ 
lution is the one where a = 2/3, for in this case the general 
solution (23) needs to be replaced by the special solution (24). 

Such a solution is shown in Fig. 2, which reports the evo¬ 
lution of e for a = 2/3 and where different lines refer to dif¬ 
ferent levels of the initial magnetization: ao = 0.5 (black 
solid line), ao = 1 (dotted light-blue line), and ao = 2.0 (red 
dashed line). As already anticipated in the previous section, 
for t > To the log term always reduces the value of e, leading 
to a faster decrease of the energy density when compared with 
the ideal-MHD limit (this is shown with a green dot-dashed 
line in Fig. 2). Furthermore, it is also clear that larger values 
of a o will lead to a faster decrease in e, as shown in Fig. 2 
(note that for ao = 2, e = 0 at r ~ 2.7 fm). 

Finally, we consider in Fig. 3 the evolution of normalized 
fluid energy density e in the case a = 2. Also in this case, 
different lines refer to different levels of the initial magnetiza¬ 
tion, ao = 0.01 (black solid line), ao = 1 (light-blue dotted 
line), and ao = 10 (red dashed line). Because the second 
term in Eq. (23) is always positive, the evolution of the energy 
density is expected to be slower than in standard Bjorken flow 



FIG. 2: Evolution of the normalized fluid energy density e for 
a = 2/3. Different lines refer to different levels of the initial mag¬ 
netization: ao = 0.5 (black solid line), cro = 1 (dotted light-blue 
line), and a o = 2.0 (red dashed line). Note that the decrease of the 
energy density is always faster than in the ideal-MHD case (green 
dot-dashed curve). 



r[fm] 

FIG. 3: Evolution of normalized fluid energy density e for a mag¬ 
netic field with a power-law decay a = 2. Also in this case, differ¬ 
ent lines refer to different levels of the initial magnetization, ranging 
from ao = 0.01 (black solid line), ao = 1 (light-blue dotted line), 
and ao = 10 (red dashed line). Note the initial “heat-up” in the case 
of large magnetizations. 


(see also Fig. 4). At the same time, the second term in Eq. (23) 
is not a monotonically decreasing function of r. As a result, 
it may produce even a temporary increase in the fluid energy 
density. This increase, which can be associated with a resis¬ 
tive “heating up” of the fluid, will depend on the values of ao 
and a and will be larger for larger values of the latter. This is 
clearly shown in Fig. 3 for r < 1 fm in the case cr 0 = 10.0; 
after this time the evolution of the energy density is monoton- 







6 



T[fm] 

FIG. 4: Evolution of the normalized energy density e in the different 
cases and when the initial magnetization is set to cro = 1. Different 
lines refer to the evolution for a = 2/3 (black solid line), a = 1 
(light-blue dotted line), and a = 2 (red dashed line). Clearly, e 
decreases more rapidly for a = 2/3 when compared to the case 
a = 1 , whereas for a = 2 it initially decays more slowly and then 
decays asymptotically at the same rate as for the ideal-MHD case 
a = 1. 


ically decreasing and asymptotically dominated by the term 

~ - 4 / 3 . 

As a way to summarize the various results presented so far 
we show in Fig. 4 the evolution of the normalized energy den¬ 
sity e in the different cases but keeping the initial magnetiza¬ 
tion fixed to cro = 1. More specifically, we show the evolu¬ 
tion for a = 2/3 (black solid line), a = 1 (light-blue dot¬ 
ted line), and a = 2 (red dashed line). Clearly, e decreases 
more rapidly for a = 2/3 when compared to the case a = 1, 
whereas for a = 2 it initially decreases more slowly and then 
decays asymptotically at the same rate as for the ideal-MHD 
a = 1 case. 


V. CONCLUSIONS 

Driven by the interest in exploring the effects of strong 
magnetic fields in the hydrodynamical description of relativis¬ 
tic heavy-ion collisions, we have studied the evolution of the 
fluid energy density following the instant of the collision and 
considering an ultrarelativistic fluid with EOS p = c 2 s e. Be¬ 
cause we are mainly interested in finding analytic solutions, 
our setup is somewhat idealized and we have therefore con¬ 
sidered one-dimensional, longitudinally boost-invariant flow 
with transverse magnetic field, i.e., a transverse-MHD flow. 
When no magnetic fields are present, this flow is known as 
the Bjorken flow [26] and although it represents a simplified 
prescription, it has served to gain significant insight on the 
dynamics of nucleus-nucleus collisions. 

We have first considered the dynamics of a one-dimensional 
MHD flow in the limit of infinite electrical conductivity and 


found a somewhat surprising result, namely that in the ideal- 
MHD limit the Bjorken flow applies unmodified. The evo¬ 
lution of the fluid energy density, in fact, is regulated by the 
same equation found for Bjorken flow and thus with an ana¬ 
lytic decay in proper time as r -4 / 3 . Of course large values 
of the initial magnetization will change the values of the total 
energy density, but the evolution of the fluid energy density 
will not be modified because of the passive role played by the 
magnetic field in this regime. This result widens considerably 
the range of applicability of the Bjorken model and shows that 
it can be used, unmodified, also to describe collisions in trans¬ 
verse MHD. 

We have also considered the cases in which the magnetic- 
field evolution is not the one prescribed by the ideal-MHD 
limit but, rather, follows a power-law behaviour in proper time 
with exponent a. The solutions in this case need to be distin¬ 
guished between the scenario in which the magnetic field de¬ 
cays more slowly than in the ideal-MHD case, i.e., for a < 1 
and when the decay is more rapid, i.e., for a > 1. In the 
first scenario, which could be realized when remnants of the 
colliding nuclei slow down the decay of the magnetic field, 
the decay of the energy density is faster. Furthermore, the 
rate at which this decay takes place is determined entirely by 
the level of the initial magnetization and modeled in terms of 
the dimensionless magnetic-to-fluid energy do = //e () . In 
the second scenario, which could be associated to a resistive 
regime in which magnetic field energy is converted to fluid 
energy via resistive losses, the evolution of the energy den¬ 
sity is more complex. In the initial stages of the evolution, 
in fact, the fluid energy density may increase as it would in 
terms of a resistive “heating up” of the fluid. The amount of 
this increase depends on the magnetic field strength and dis¬ 
sipation and hence will increase with op and a. However, as 
the fluid further expands, its energy density will decrease with 
an asymptotic rate that is the same as in the Bjorken flow, 
i.e., oc t- 4 / 3 . 

The work presented here could be extended in a number of 
ways. First, one could search for analytic solutions in one¬ 
dimensional MHD in a Landau-type flow scenario. Second, 
one can consider an explicitly finite electrical conductivity 
as the simplest model for a one-dimensional MHD flow with 
chiral fermions. Results on these topics will be presented in 
forthcoming papers. 
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Appendix A: Covariant derivative for Bjorken expansion 

In this appendix we sketch briefly the steps that are needed 
to derive the energy-and momentum-conservation equations 
discussed in the main text. Before we start, we should 
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comment about the notation used in defining the energy- 
momentum tensor (1). Bearing in mind that in general- 
relativistic calculations the standard choice for the signature 
is a spacelike one, i.e., (—, +, +, +), the energy-momentum 
tensor in GRMHD is normally defined as [31, 32] 


= (e + p + b 2 ) // + 



_ b n b » s (Al) 


with 7 = 1 /\/l — v 2 the the Lorentz factor and B is the mag¬ 
netic field three-vector measured by an Eulerian (or normal) 
observer. The modulus of b 1 ' is then given by 


b 2 = b^ = 




(A3) 


where b is the magnetic field four-vector measured in a co¬ 
moving frame and has components given by 


(7 v ■ B, -1- 7 v ■ B v 


(A2) 


With this clarification in mind, we go back to our special- 
relativistic setting with energy-momentum (1) and consider 
the projection of the energy-momentum conservation equa¬ 
tion dJT = 0 along the fluid four-velocity u which reads 




u il u^u v d v (e + p + B 2 ) + (e + p + B 2 ) u^d v - u^d u 


B- 

p + -j-1 r 


o 

— Utl d v {B^B U ) = 0 , 


i v d„ (e + p + B 2 ) + (e + j) + B 2 ) d v u v - u v d v ^p + + B^B v d u u^ 

a,(e + P + s») + t±i!±^!-a,( P+ 2!) 


r>2 

dr[e+ — 


e + p+ B 2 


(A4) 


where we have used Eq. (1), u^B^ = 0 and = 0, Similarly, the projection of the conservation equation 

since = (uq,0, 0,w 2 ) and B M = (0, B x , B y , 0) in our d^T^ 1 ' = 0 in the direction orthogonal to the fluid four- 


transverse-MHD setup. 


velocity gives 


K v d a r av = o, 

(e + p +B 2 )h IJI yd a (u a u’') - h llv d v ^p + - h^ v d a (B a B’') = 0, 

(e +p + B 2 )u a d a u,j, - h /tM d ,J ^p + - B a d a B fl - B^d a B a + ii^d a (B a B v ) = 0, 

(e + p + B 2 )d T u l _ l - h ll „d 1 ' fp + ) ~ B a d a B fl - B,j,d a B a - u^B a B , 'd a Uv = 0, 


(A5) 


where we have introduced h as the orthogonal projector to u, 
i.e., h u = 0, where h^ a = rj„ v — u^u u . The last three terms 
vanish, because B M is assumed to be constant in transverse 
direction. This then leads to Eq. (7). 


Appendix B: Frozen-flux theorem 

In this appendix we show that the evolution of the magnetic 
field and of the entropy density are strictly related in the ideal- 
MHD limit. The arguments reported below are well known 
and can be found in a number of textbooks (e.g.. Refs. [33, 
35]), but we recall them here for completeness. We start from 
the definition of the covariant (or Lagrangian or convective) 
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time derivative given by 


D _ d 
Dt dt 


+ u ■ V, 


(Bl) 


where u is the fluid velocity. If x is the position of a fluid 
element, this will be advected with the flow and hence have 


Dx 

lot 


= 0 . 


(B2) 


However, if £ is a vector separating two fluid elements at a 
given instant, the corresponding Lagrangian derivative will 
not be necessarily be zero, but is actually given by 


^=lVu. (B3) 

Stated differently D£/Dt = 0 only for a fluid with uniform 
velocity u. In all other cases, the vector £ will change its 
length and/or orientation in the presence of a velocity gradi¬ 
ent. 

Next, we consider the conservation of rest mass, which can 
be expressed as 


^ = -pV • u, (B4) 

where p is the rest-mass density of the fluid. In ideal MHD 
the induction equation takes the well-known form 


dB 

~dt 


= V x 


[u x B'j 


(B5) 


and the frozen-flux theorem states that the magnetic field lines 
are frozen in the fluid and can be identified with the worldlines 
of fluid elements. To see this, we use the following vector 
identity 


Vx^ux B j = B-Vu—B • uj -u-X7B+u • Bj , 

(B6) 

together with V ■ B = 0 in Eq. (B5) to obtain 

^ = B • Vu - B (V -uj . (B7) 


Together with the conservation of mass, the above equation 
can then be written as 


D 

Dt 




(B8) 


This is exactly the same equation satisfied by the separat¬ 
ing vector £ [Eq. (B3)]. Therefore a magnetic field line is 
advected and distorted by the fluid in the same way as a fluid 
element. If the fluid expansion takes place isentropically, the 
total entropy of the system remains constant and the entropy 
density s will satisfy the same conservation equation as the 
rest-mass density, i.e., 


/ J s —* 

— = -sV • u . (B9) 

Dt 

From the arguments made above, it follows that the quantity 
B/s will behave as the quantity B/p and hence satisfy the 
equation 


D 

Dt 




(BIO) 


which is identical with Eq. (B8) except p is replaced by s, 
i.e., for this case we also have the magnetic fluxes frozen in 
the system. 


Appendix C: Limit for the log term 


In this appendix we discuss how to evaluate the second term 
in Eq. (23) in the limit in which a —> (1 + c^)/2, i.e., the limit 


lim 

a—H l+c|)/2 


i T o/ T ) +Cb ~ fa fry 
1 + c? - 2 a 


(Cl) 


We first increase the power exponent a by an infinitesimal 
amount e > 0 and then take the limit e —► 0. In this case, 
Eq. (Cl) becomes 


, im , 

e->o 1 + Cg — 2 (a + e) 


(C2) 


and setting a = (1 + c 1 2 3 s )/2 we obtain the desired result 



lim 

e ->0 


1 - (W t ) 2£ 

(—2e) 
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